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1. Introduction 

Nonperturbative studies of QCD at finite baryon density and finite temperature is an active 
field of investigation on the lattice with implications for current experimental studies of the 
transition between the hadronic phase and quark-gluon plasma phase [1, 2, 3, 4]. Due to the 
well known "sign-problem" associated with lattice simulations at finite chemical potential, 
a large fraction of the numerical simulations compute coefficients in the Taylor expansion in 
chemical potential [5, 6]. Typically one computes susceptibilities (second order coefficients) 
but higher moments have also been computed [3, 7]. The Taylor expansion is expected to 
diverge close to the critical point where fluctuations dominate the physical behavior [8]. 

We will focus on the quark number susceptibility for definiteness in this paper but the 
formalism in this paper carries over to all coefficients in a Taylor expansion. Let Z(fi, Nt) 
be the grand canonical partition function on a iV 3 x Nt lattice at a given chemical potential 
/i. The dependence of the partition function on other quantities like gauge coupling, quark 
masses etc. have been suppressed. We also assume that we will take the thermodynamic 
limit, N — > oo, and 1 will play the role of temperature on the lattice. The quark number 
susceptibility is defined as 



,. r , .. 1 d 2 In Z(n,N T ) 
X{N T ) = lim 



JV^oo iV3jV T dp? 



(1.1) 

/x=0 
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This quantity has been studied both on the lattice [1, 9] and in other non-perturbative 
approaches [10]. 

A naive introduction of the number density operator on the lattice leads to divergences 
at zero temperature in the following sense. A proper definition should yield a result that 
goes as j^j- and the presence of a non-zero x(oo) implies a divergence in the continuum 

limit. Such a divergence can be removed by replacing U±(x) and U\{x) on the lattice by 
e^U^x) and e~^Jj\{x) where U±{x) is the link variable on the lattice connecting x and 
x + A [11] and will be referred as the Hasenfratz-Karsch (H-K) prescription. Let T4 denote 
operator 

T^(x) = U 4 {x)ip{x + 4). (1.2) 
For naive fermions, the source term added to the action would be 

{^74^ for naive insertion; 

"074 — — 1 ^ >Ti ^ e 1 ^ L/ ifj for H-K insertion. 

The naive insertion would result in a non-zero x(°°) but the in H-K insertion it would lead 
off as -^2-. 

The aim of this paper is to address the introduction of the chemical potential into the 
overlap Dirac operator [12]. It was introduced using the H-K prescription by Bloch and 
Wettig [13]. The Wilson-Dirac operator, H w (fi), that appears as the kernel in the definition 
of the overlap Dirac operator is not Hermitian in the presence of a chemical potential. The 
massless overlap Dirac operator in the presence of the chemical potential as per the H-K 
prescription is 

DbwM = 1 + 75£ f- M) - (1-4) 
Assuming H w {fj) is diagonalizable in the form 

H w (fj) = RAJRT 1 , (1.5) 

with A being a diagonal matrix with complex entries, one sets 

e(H w (fi)) = Re(ReA)R-\ (1.6) 

Such a choice can be justified using the domain wall formalism [14] since e(H w ) is the 
limit of tanh(L s H w ) with L s being the extent in the extra continuous dimension [15]. In 
spite of the fact that V(fi) = ^^e{H w {fi)) is not a unitary operator for fi^O, the massless 
propagator defined in the usual manner [15, 16] as 

is still chiral: Let us write 

*-(" ]); A.( A + M; ReA ± >„. (1 .8) 
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Then one can show that 



i 



G(/i) 



( 





j3a 



-i 







) 



(1.9) 



showing the chiral nature of the massless propagator assuming that a and 5 are invertible 
complex matrices. Both the energy density [17] and the quark number susceptibility [18] 
for the case of free overlap fermions were analyzed within the Bloch-Wettig scheme and 
shown to have the correct continuum behavior. One possible objection to the Bloch-Wettig 
scheme would be that the chemical potential is coupled to all the fermions from the view- 
point of domain wall fermions and not just to the chiral fermions. This objection could 
also be raised for Wilson fermions where the chemical potential is coupled to the massless 
mode as well as the massive doublers but we should keep in mind that we have infinite 
number of fermions in the case of domain wall fermions. Since we have an explicit operator 
for the physical fermions, it should be possible to couple the chemical potential just to 
the physical fermions in the overlap formalism. Another possible objection already raised 
in [13] is that the number of topological modes can depend on the chemical potential since 
the topological charge is 



where tr is the trace over the spin components at a fixed point x. Since e(H w (fi)) depends 
on //, Q can change with [i in the Bloch-Wettig scheme. In fact q(x) was shown to depend 
on \x in [19]. 

We will start from the first principles of the overlap formalism [20] and add a chemical 
potential by coupling it to the number density operator. We will consider the general 
case where the chemical potential for the fermions with positive and negative chirality 
can be different. We will make a class of choices for the number density operator all of 
which reduce to the standard continuum limit. We will focus on the case of free overlap 
Dirac fermions and show that divergences in the quark number susceptibility is a generic 
feature. Since every member in our class of choices will contain the correct in the free 
quark number susceptibility, we will be able to write down a large class of number density 
operators that will not contain a divergent piece. 

As a first step, we will review the generating functional 1 for a massless vector like gauge 
theory. None of this is new and can be found in one form or other in [15, 20]. In particular, 
we will explicitly write down the chiral transformations in the overlap formalism and derive 
the massless overlap Dirac operator [12]. This will help us introduce the chemical potential 
and derive the associated generating functional for the massless overlap Dirac operator 
with a chemical potential. The resulting class of operators will be different from L>bw(a*) 
in (1.4) but the generating functional will be invariant under chiral transformations. All 
these operators will have two features that separate them from the operator in the Bloch- 
Wettig scheme: Chemical potential will only be coupled to the physical chiral fermions and 
topological charge will not depend on the chemical potential. Finally, we will discuss the 

1 We will only consider the generating functional for a fixed gauge field background through out this 
paper. 




(1.10) 
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addition of a mass term to the generating functional that destroys the chiral symmetry. 
We will end the paper with a discussion of the quark number susceptibility for free overlap 
Dirac fermions arising from our generating functional. 

2. The overlap formalism 

This section is essentially a repeat of what can be found in [20] . The main difference is that 
we only use one Hermitian Wilson-Dirac operator in the construction of overlap fermions. 
This amounts to setting the Wilson mass parameter to infinity on the inactive side. This 
changes the details as it pertains to the generating functional and since we will derive the 
generating functional for massive overlap fermions in the presence of the chemical potential, 
we will need explicit results from this section. 

At the heart of the overlap formalism is the massive Hermitian Dirac operator which is 
usually realized on the lattice using the Wilson-Dirac operator, H w . The notational details 
are given in Appendix A. The generating functional for a vector like theory with massless 
fermions is given by [20] 

Z&0 = R(-\^ RdR+ t RU «\+) R L (+\^ Ld ^ LUL \-)L. (2.1) 

Dirac valued operators, a R) L,a R L , obey canonical anti-commutation relations separately 
for the R and L sets. We write, 

«»-(%)• (2 ' 2) 

in terms of their respective Weyl components. The sources, 

£ = ex); Z=(f L ) ( 2 - 3 ) 

are Dirac values Grassmann variables that couple directly to the physically relevant fields 
and they anti-commute with the fermionic operators. \+)R t L are lowest states of 

Hr,L = -a RL H w a R ^ L (2.4) 
and |— )r j l are lowest states of two many body operators, 2 

Tr,l = -a RtL 75aR,L- (2.5) 
We now make several remarks on the generating functional. 

1. The phase choice for are tied together since they are the ground states of 

identical many body operators. The same is true for | — Therefore, the generating 
functional is unambiguous and does not not depend upon the phase choice present 
in the unitary matrix, U, that diagonalizes H w . 

2 The pair Hr and Tr replaces the H± in [20]. That one can take the Wilson mass parameter to infinity 
on one side follows from the discussion of the topological charge in section 8 of [20]. 
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2. It does not depend on the ordering of the operators since the two terms in the 
exponent commute with each other in both factors. 

3. It is clear that cLr and u R are the propagating degrees of freedom in the first factor 
since r(—\ur and r(— \d R are both zero. The converse holds for the second factor. 

4. The generating functional is invariant under global chiral transformations: 

tR^e^te & "> Cfle-**; £ L W^£ L ; h^he~ itpL . (2.6) 



5. The explicit expression for the generating functional can be obtained by small mod- 
ifications to the details presented in [20] and the result is 



a 



det or 



In particular, the propagators for right-handed and left-handed fermions are 
R(-\u Rj d Ri \+) R 



Gl 



r(-\+)r 

and they obey the relation 



ij ' L 



h{+\d ] Lj u Ll \-) L 



■A+\-)i 



a 



(2.7) 



(2.8) 



(2.9) 



In practice, one can avoid exact diagonalization of H w which is needed for the compu- 
tation of U and (2.7) since one has an overlap-Dirac operator for vector like theories [12]. 
Consider the unitary operator, 

V = l5 e(H w ). (2.10) 

It follows from (A. 6) and (A. 7) that 



l + V 



u 



a 



l-V, 



and therefore 



G Q 



5 J ' 2 
-V ( 



-U 



7 




l + V \f3a~ 1 
is the massless overlap Dirac propagator. Since [20] 

det a 



detU 



we have the identity, 



and therefore, 



det <yt ' 

det a det ofi = det 5 det 6' , 
l + V 



det 



det a det or . 



The massless overlap Dirac operator is given by 



(2.11) 
(2.12) 

(2.13) 
(2.14) 
(2.15) 

(2.16) 
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Our generating functional in (2.7) can be written as 3 

Z(£ = det D e- f ^ G ^. (2.17) 

We can use (2.12) and (2.15) along with an efficient implementation of V to compute the 
generating functional. Note that the operator in (2.16) is not identical to the operator 
used to compute the propagator in (2.12). This is an essential ingredient of the overlap 
formalism. 



3. Introduction of the chemical potential 

Consider the generating functional 

z&b m,M = H (-| e ^ d «+^4+«l i A«'ii i | + ) ili ( + | e ^4+^^-4A i ^|_) i) (3.!) 

where (ir(iir) and /ti(/Xi) are operators that parametrically depend on the chemical po- 
tentials, fin and fiL- 

In order to obtain a formula for the generating functional, we start by noting that 

f d(Rd C Re -^-M, = A^ (3.2) 

and 

/ d( L d( L e-^ u ^ d i~^ = e -4Ai«£ j (3.3) 

where (r,Cl,Cr> Cl are Grassmann variables that anticommute with all fermionic operators 
and Grassmann variables. Therefore (2.1) and (2.7) gets modified to 



Z fiR, fiL) = det a deter 

dC L dCLe-~ CLa+ ^ L -^ L ^ a ^^ L - CL) 
det a det (l + /i/j/3a _1 ) det a 1 det ( 1 — fii, [/3~~ 



a/J-l+AH M" 1 ] 1 -*L ^ . (3.4) 

Note that the generating functional is invariant under the chiral transformation given 
in (2.6) since the introduction of the chemical potential does not mix the two chiral sectors. 
To be consistent with the continuum definition of the chemical potentials, we require 

Afi(0)=0; A-l(aO = -Ajr(-aO- (3-5) 

Note that the fermion determinant is not real and positive for a real quark chemical 
potential, fin = fi^ = fi, but it is real and positive for an isospin chemical potential, fiR = 
—fiL = /-*■ These are standard properties in the continuum that are correctly reproduced 
by the overlap formalism. 



3 The presence of —174 is due to our choice of basis for the chiral sources. 
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There are several options for Ar(aO- We will address this issue when we analyze free 
fermions. 

We can write (3.4) in a compact form by introducing the massless overlap Dirac oper- 
ator and propagator in the presence of a chemical potential. If we define 

flR 

,H 



N 



'1 + V 



+ N 



1-V 

then a straight forward computation shows that 

Z(Z,£;p,R,fi L ) = det£ (£ fl ,£z)e-^ G °^^. 



(3.6) 
(3.7) 



4. Introduction of the fermion mass 

A mass term can be added in analogy with the chemical potential term [20] . The generating 
functional in the presence of a chemical potential and mass term is defined to be 

Z(€,&frR,P>L]Tn) = [r(-\®l{+\] 

e iRd R +i R ul i +ul l il R d R +^ L dl j +^ L u L ~d\ j il L u L +md^d R -m,ul i u L 

[\-)l®\+)r] (4.1) 
One can show using methods similar to one in (3) (see Appendix B for details) that 

Z(£,£;fiR,fiL]rn) 
= det adet (l + (LRfioT 1 ) det det (l — jli (/3a -1 ^ 



det 



1 + mf 



1 



a/3 1 + pi R 



1 



(a/3- 



-nt 



m2 + ( a /3-l+A fl )((c« / 9- 1 ) t -ML) 



5B+a(a / 3- 1 +Afl' 



m 2 +((a5- 1 ) t -Ai)(«(3 _1 +Ail) 



2 + (( Q /5-l)t- Ai )( a/ 3-l +Ai j) 



2 + (^-i+A ii )((«/3- 1 ) t -A L ) 



(4.2) 



We can write (4.2) in a compact form by extending the massless overlap Dirac operator 
in the presence of a chemical potential define in (3.6) to include a mass as follows: If we 
define 

ATf \ ( 77 

N[m) ■ 



(4.3) 



fjL L m 

1 + V 1-V 

D (fl R ,p, L ;m) = — \-N(m) — - — 



2 

1 + V 



1-V 



+ N(m) 



G (fiR, fiLirn) 
then a straight forward computation show that 

z{l Ah, Al; m) = det z> (Afl, £l; m) e -** G ^*^ m K 



(4.4) 



Note that the expressions reduce to the usual ones [16] for the massive overlap Dirac 
operator without a chemical potential. 
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5. Free overlap fermions 



In order to work out the energy density and quark number susceptibility for free massless 
quarks, it is best to work in momentum space. Let us assume that we have converted 
to creation and annihilation operators in momentum space by the appropriate unitary 
transformation. We will work on a oo 3 x Nt lattice. The allowed spatial momenta are 
Pk € [— vr, 7r] . The Matsubara frequencies in the Nt direction are oj n = ± 2 ™^~ 7r ; n = 
0, • • • , ^ - 1 and we will assume that Nt is even. The hermitian Wilson-Dirac operator 
for a fixed momentum takes the form 

B. = fc c ) I") 

\c l —w J 

where 

w = 1 y sin — + 2 sin — m w ; c = i y smpk — sino; n ; 

k k 

= = sin 2 pk + sin 2 uj n . (5.2) 



The positive eigenvalues come in one doubly degenerate pair, A = Vw 2 + cct , per momen- 
tum block and the corresponding pair of orthonormal eigenvectors are 

, 1 (V V (5.3) 

Note that A and w are even functions of oj n at a fixed pk- Assuming a vector like chemical 
potential, we write 

\-n*(-n,u n ) J 
per momentum block where n(0,uj n ) = 0. We write, 

n(fjb,LJ n ) = n r (fjt,u} n ) + ini(fi,u n ) (5.5) 

where n rj i(fi,u n ) are both real functions. The free energy density is given by 

3 

]nZ(»,N T ) = f[[ ^l^{_21n(2A) + ln( 5 2 + / l )} (5.6) 

where 

g = X + w + n((jL,u n )n*(-fj,,u n )(\ - w) - [n(p,uj n ) + n*(-fi,u n )] smuj n , (5.7) 

and 

3 

h= [n(n,u n ) - n*(-fi, w„)f ^sin 2 p fc . (5.8) 

k=i 
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The quark number susceptibility is 

1 fx dpi s-^ 4 (^ 2 ~ r2 ) ( sin2 Wn ~ Sfc=i sin2 Pfc) - 4s(w n )(A + w) sine 



x(N T ,m v 



where 



q(oj n ) 



i=l 



n/ t;t. 



[X + wf 



dni{n,u n 



r(cj n ) 



dn r (n,w n ) 



s(w„) 



<9 2 n r (/i,u; n ) 



d 2 fi 



(5.9) 
(5.10) 



/i=0 



The sum over n can be non-zero only if 4(g 2 — r 2 ) is an even function of Lo n and s is an 
odd function of u n . 

Two choices one typically makes are 



2m„ 



for naive insertion; 



(e^-DrHe-^ for H-K insertion; 



(5.11) 



resulting in 



n(ji,u r , 



1 f i/i for naive insertion; 

2772^, \ sincj n — sin(w n — i/j,) for H-K insertion. 



(5.12) 



The -7— factor comes from the tree- level wavefunction renormalization [161. For the naive 
insertion, 

L - ; r = 0; s = 0. (5.13) 



2m„ 



For the H-K insertion, 



COS UJ r 

2m,,, 



r = 0; s 



sin uj Ti 
2m,„ 



(5.14) 



We can make a modification to the H-K insertion such that s = but keep q to be the 
same. This corresponds to 



Ak(m) = 7^— sinh |(T 4 + t\) => n(/x, u;„ 



Keeping only the forward derivative will result in 



1 Ju n r ---l_ M 



— cos uj n smn — . 



with 



n(n,u n ) = —e luln sinh 
m,„ 2 



cos w n sm w n 
q=— ; r = — ; s = 0. 



This leads us to consider the restricted class of operators, 



Ar(^) = ZT~ sirm T7 T 4 ^ nJ Wri 



■;??,, 



— e ywB sinh^, 
m w 2 



(5.15) 

(5.16) 
(5.17) 

(5.18) 
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with 



= cosjujn r= sinjo; n _ s _ q. g 2 r 2 = cos 2 3^n (5 19) 

2m w ' 2m^ ' ' (2m w ) 2 



The operator T\ is the j-th power of the operator T4 defined in (1.2). This motivates us to 
study 



1 JL /-t cos(jw n ) (^sin 2 w n - Ylk=i s[n Pk 



i=l"-« n 



(X + w) 



-, j = even. 

(5.20) 



and we conjecture that 



Ij(N T) m w ) = Ij{oo, m w ) + -^ + (Nf 4 ) . (5.21) 

We are unable to prove (5.21) analytically since we need to work at a finite but large 
Nt- The resulting contour integral is complicated due to the presence of the branch cuts 
associated with square roots signs in definition of A. We resort to a numerical check of the 
conjecture. For this purpose, we convert the three dimensional integral over the momenta 
to a three dimensional sum with N points in each direction. Let Ij(N, NT,m w ) denote 
this sum. Keeping N fixed, we compute the Ij(N, 00, m w ) by computing the quantity for 
a very large Nt- We take this result and extrapolate to N — > 00 and obtain Ij(oo,m w ). 
We then consider Ij((Nt, NT,m w ) — Ij(cc,m w ) with £ = 3, 4 in order to extract the 
coefficient in (5.21). We numerically evaluated 



T 



RjiCNnmv) = 3N T (I,{CNT,N T ,m w ) - ^(00, m w )) - 1 (5.22) 

for j = 2,4 at £ = 3,4 and m w = 1,1.5. It should approach zero with corrections of the 
order of Fig. 1 shows numerical evidence in support of the conjecture: The fits show 

that Rj(4,Nt, 1) approach zero as but the corrections get larger with j. Fig. 2 shows 
that Q = 4 is sufficiently large for the numerical computation at hand and Fig. 3 shows 
that the effect of the Wilson mass, after taking into account the tree level wavefunction 
renormalization factor, is minimal. 



6. Conclusions 



We have discussed the problem of introducing the chemical potential within the overlap 
formalism. The formalism presented here couples the chemical potential only to the phys- 
ical chiral fermions and the overlap definition of the topological charge is unaffected by 
the value of the chemical potential on the lattice. We have discussed a large class of op- 
erators coupled to the chemical potential. Each one of them results in a divergent quark 
number susceptibility but all of them have the correct continuum limit after the divergent 
contribution is subtracted. We envision using the formalism developed in this paper in two 
ways: 
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Figure 1: Numerically obtained plot in support of the conjecture. 

1. Since the main physics aim in the Taylor expansion method of studying the phase 
transition from the hadronic phase to the quark-gluon plasma phase is to look for 
growing fluctuations close to the phase transition, we could consider, say the quark 
number susceptibility, by taking the difference between its value on a finite tempera- 
ture lattice and zero temperature lattice for the same value of lattice gauge couplings 
and lattice quark masses [19]. This will remove the zero temperature divergences and 
enable a proper study of the fluctuations close to the transition temperature. The 
chemical potential was coupled to only those physical fermions that are confined to 
the 4D domain wall at the origin of the fifth dimension in [19]. However the num- 
ber density term did not commute with the overlap Hamiltonian. In this paper the 
chemical potential is coupled to the conserved number density operator explicit in 
the many body overlap formalism. 

2. Consider the linear combination, 

oo 

= 5^ c jAk(a*)- (6-1) 

3=1 

The associated free fermion quark number susceptibility for this choice of the number 
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Figure 2: Numerical evidence that £ = 4 is large enough for the computation at hand, 
density operator will be 

oo 

x{N T ,m w )= ^2 c ii c i2 / ji+j 2 ( Ar T,^)- (6.2) 

As per the conjecture in (5.21), the quark number susceptibility will give the correct 
behavior as long as 

5> = 1. (6.3) 

i 

In order to have no divergence, we need to choose the coefficients, Cj, such that 

oo 

c h c n I h+j2{°°^ m w) = 0, (6.4) 

h ,32=1 

and this is only one condition and we have infinite coefficients. Therefore, we can 
find a large class of number density operators that have the correct finite behavior in 
the free fermion limit. Addition of gauge fields cannot give rise to new divergences as 
long as all couplings and masses are properly renormalized. We could therefore choose 
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~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 
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• m =1 

w 
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Figure 3: The effect of the Wilson mass parameter on the subleading terms is minimal. 

c\ = a and C2 = (1 — a) with a chosen to cancel the free fermion divergence. This 
will provide one good choice for the number density operator in the full interacting 
theory. There is clearly one weakness to this approach since one has to tune a to 
cancel the divergence. 

Dynamical simulations of overlap fermions are currently being performed both at zero 
temperature [21] and finite temperature [22]. It would be interesting to compute the quark 
number susceptibility in such simulations using the operators presented in this paper and 
compare them with the complimentary approach of [13]. The two approaches are quite 
different in their construction. We only couple the chemical potential to the physical 
fermions. In addition to coupling the chemical potential to the physical fermions, it is 
also coupled to an infinite number of regulator fields in the Bloch-Wettig scheme. A 
related issue is the definition of the topological charge. It does not depend on the chemical 
potential in our scheme contrary to the Bloch-Wettig scheme. The continuum definition of 
the topological charge based on the counting of the zero modes of the chiral Dirac operator 
is not affected by the insertion of the chemical potential and our scheme maintains this 
continuum property on the lattice. Since we expect to work with large values of lattice 
chemical potential particularly close to the physical transition that separates hadronic 
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matter from quark-gluon plasma, the lattice spacing effects arising from the topological 
charge depending on the chemical potential could be severe. In spite of these differences, 
we expect lattice simulations to show that the two approaches agree in the continuum limit. 



A. Details of the Hermitian Wilson-Dirac operator 



Ct —B + m 

be the massive Dirac operator with < m w < 2. In the above equation 

(T^)(x) = U^x)<f>(x + ft) (A.2) 

with 

n ~(i I)'- ° 2 = {°i oO i8 »-(J "'-{l, ") (A - 3) 

in d = 4 and 

oi = 1; a 2 = -i (A.4) 
in d = 2. The Dirac matrices in the chiral basis are 

-a y). *-(; -0- 

We will assume that H w is a 2n x 2n matrix and we will assume the gauge field is in the 
zero topological sector. 4 Let 

H W U = UK, (A.6) 



with 



> T -[°p ]); A = diag(A+ •••,A+ -\-,...,-\~). (A.7) 



and A^ > for all i. 



B. Details of the derivation of generating functional for massive overlap 
fermions with a chemical potential 

Like in (3), we can introduce Grassmann variables, Cr>Cl,CriCl, that anticommute with 
all fermionic operators and Grassmann variables to rewrite (4.1) in terms of (3.1) as 

Z{li-fi R r^m) = I d( L d( R d( L d( R e^-^ 



4 The overlap formalism is designed to work in all topological sectors and this is one of the important 
features of the formalism. But, we can restrict ourselves to the zero topological sector to simplify the 
discussion. All results will trivially extend to all topological sectors. 
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\ R l—\ (g) e-V^Cfl^fl-V^C-ffUB-V^C-LML-V^CL^l 

Using (3.4) the above equation can be written as 



^RdR+^RU^+U^flndR+^Ldl+^LUL-dljlLUL n_\ L (g) . (B.l) 



.£(£,£; pLR,fi,L', Tn) = det a det (l + fiRpa r ) det det ^1 — jxL (P 

dC L dQ R dC L d(; R e^L-ZL(R 



e (gR-V^Cil) a ^-l +Afl (€fl-V^Cfl) 



and 



and 



l+m 2 GiGT l+m 2 GjGi 



and it follows that 



from which the last two lines of (4.2) follows. 



Acknowledgments 



(B.2) 



We can write the exponent in the integrand as 

S = (F( - ^ {(G£ + &() + ^ (B.3) 

where 

C = (Cr Cl); C=(^); (B.4) 

Gl = a/3-i 1 +AH ! ^ 2= (a/3-i)t-AL- ^ 
Integration of £ an d C yields 

detFe«( G - mGF " lG )«. (B.7) 
The second line in the last equality of (4.2) is det F. One can show that 

mGt 1 — 



F -l = I 2 l+m 2 G!G^ l+m 2 G'Gi j ^ 



iWG.Gt 1 l+m 2 G|Gi 2 j (B . 9) 

2 l+m 2 G 1 Gl 1 l+m 2 GlGi 2 / 
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